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Abstract. The glass transition is considered within two toys models, a mean field 
spin glass and a directed polymer in a correlated random potential. 

In the spin glass model there occurs a dynamical transition, where the the sys- 
tem condenses in a state of lower entropy. The extensive entropy loss, called com- 
plexity or information entropy, is calculated by analysis of the metastable (TAP) 
states. This yields a well behaved thermodynamics of the dynamical transition. The 
multitude of glassy states also implies an extensive difference between the internal 
energy fluctuations and the specific heat. 

In the directed polymer problem there occurs a thermodynamic phase transition 
in non-extensive terms of the free energy. At low temperature the polymer condenses 
in a set of highly degenerate metastable states. 



1 Introduction 

The structural glass transition is said to occur at the temperature T g where 
the viscosity equals 10 14 Poise. The question why this transition occurs is 
often "answered" (more correctly: avoided) by saying that it is a dynamical 
transition. Surely, there is a continuum of time scales ranging from picosec- 
onds to many years; at experimental time scales there is no equilibrium. 
Nevertheless, since some 20 decades in time are spanned, one would hope 
that equilibrium statistical mechanics can be applied in some modified way. 

Crudely speaking, the observation time will set a scale. Processes with 
shorter timescales can be considered in equilibrium; processes with longer 
timescales are essentially frozen, as if they were random. To provide a (non- 
equilibrium) thermodynamic explanation of a model glassy transition will be 
the first subject of the present work. 

Intuitively we expect that the resulting free energy is given by the loga- 
rithm of the partition sum, provided it has been restricted to those states that 
can be reached dynamically in the timespan considered. This non-equilibrium 
free energy will then differ from the standard case, and need not be a ther- 
modynamic potential that determines the internal energy and entropy by its 
derivatives. 
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Experimentally one often determines the entropy S exp , and thus the free 
energy F exp = U — TS exp , from the specific heat data by integrating C/T 
from a reference temperature in the liquid phase down to T. As long as the 
cooling rate is finite there remains at zero temperature a residual entropy. In 
the limit of adiabatically slow cooling it vanishes. 

Alternatively, a glass can be seen as a disordered solid. In this description 
the liquid undergoes a transition to a glass state with extensively smaller 
entropy. These states are sometimes called "states", "metastable states", 
"components" or, in spin glass theory, "TAP-states" . As the free energy 
then becomes much larger, it is not so evident from thermodynamic con- 
siderations why the system can get captured in such a state with much and 
much smaller Gibbs- weight ~ cxp(— volume). The explanation is that the 
condensed system then has lost part of its entropy, namely the entropy of 
selecting one out of the many equivalent states. This part, J, is called the 
configurational entropy, complexity or information entropy [1] [2]. Its origin 
can be understood as follows. When the Gibbs free energy F a of the rele- 
vant state a has a large degeneracy A/a = exp(Z a ), the partition sum yields 
z = E a «xp(-/3F a ) ss M a cxp(-/3F a ), so F = F a - TX S is the full free energy 
of the system. The entropy loss arises when the system chooses the state to 
condense into, since from then on only that single state is observed. [3] As 
the total entropy S = S a + In is continuous, so is the total free energy. For 
an adiabatic cooling experiment Jackie has assumed that the weights p a of 
the states a are fixed at the transitional] which implies that the free energy 
difference between the condensed phase and the liquid is positive and grows 
quadratically below T c . This explains the well known discontinuity in quan- 
tities such as the specific heat. However, it may seem unsatisfactory that this 
higher free energy branch describes the physical state. 

We shall first investigate these questions for the dynamical transition of 
a mean field spin glass model, and then for the static transition of a directed 
polymer in a correlated potential. 



2 The p-spin glass 

We first analyze these thermodynamic questions within a relatively well un- 
derstood spin glass model, the mean field p-spin interaction spin glass. For a 
system with N spins we consider the Hamiltonian 

"H = — Ji 1 i 2 ---i p Si 1 Si 2 ■ ■ ■ Si p (1) 

il<i2<---<i p 

with independent Gaussian random couplings, that have average zero and 
variance J 2 p\/2N P ~ 1 . 

Kirkpatrick and Thirumalai[4] pointed out for the case of Ising spins that 
for p = 3 there is a close analogy with models for the structural glass transi- 
tion, and that its properties are quite insensitive of the value of p as long as 
p > 2. 



Complexity as the driving force for glassy transitions 3 

The spherical limit of this model, where the spins are real valued but 
subject to the spherical condition ^ - Sf = N, is very instructive. It has 
received quite some attention recently. The static problem was solved by 
Crisanti and Sommers. [5] There occurs a static first order transition to a 
state with one step replica symmetry breaking (1RSB) at a temperature T g . 

The dynamics of this spherical model was studied by Crisanti, Horner 
and Sommers (CHS) [6] and Cugliandolo and Kurchan [7]. Both groups find 
a sharp dynamical transition at a temperature T c > T g , which can be inter- 
preted on a quasi-static level as a 1RSB transition. This dynamical transition 
is sharp since in mean field the metastable states have infinite lifetime. For 
T <T C one of the fluctuation modes is massless ( "marginal" ) , not unexpected 
for a glassy state. At T~ there is a lower specific heat. 

These dynamical approaches are the equivalent for the spin glass of the 
mode coupling equations for the liquid-glass transitions. [8] At a critical 
temperature T c > T g a dynamic phase transition has been reported. The 
presence of a sharp transition has been questioned, however. [9] 

CHS integrate C/T to define the "experimental" entropy S exp and the re- 
sulting free energy F exp = U — TS exp exceeds the paramagnetic one quadrati- 
cally. The interpretation of metastable states ( "TAP-states" ) in this system is 
discussed in rcf. [10] The statistics of those states was considered by Crisanti 
and Sommers (CS). [11] Assuming that the result of long time dynamics fol- 
lows through being stuck in the metastable state of highest complexity, they 
reproduced the "experimental" free energy obtained of CHS. This confirms 
Jackie's prediction of a quadratically higher free energy in the glassy state. 

The long-time dynamics of 1RSB transitions fixes qo (= in zero field), q\ 
and x, which are just the plateau values and the breakpoint of a related Parisi 
order parameter function, respectively. In p-spin models they can simply be 
derived from a 1RSB replica calculation provided one fixes a: by a marginality 
criterion for fluctuations on the q± plateau. [4] [6] [7] The present author 
recently assumed that this is a very general phenomenon. [12] This was 
motivated by the expectation that a dynamical transition will automatically 
get trapped in a state with diverging time scale, if present. In a Potts model 
this then predicts a dynamical transition with marginal qo plateau and stable 
qi plateau. 

As the replica free energy is minimized in this procedure, it lies near T~ 
below the paramagnetic value and has a larger slope. Though this is exactly 
what one expects at a first order phase transition, it is a new result for 
dynamical glassy transitions. It is the purpose of the present work to discuss 
the physical meaning of the mentioned free energies. 
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2.1 The replica free energy 

At zero field the 1RSB replica calculation involves the plateau value qi and 
the breakpoint x. It yields the free energy [5] 



AT 4 4 

-^log(l-^i) + glog(l-gi) 

where £ = 1 — x. The first term describes the paramagnetic free energy. Here 
and in the sequel, we omit the T = oo entropy. It is a constant, only fixed after 
quantizing the spherical model, [13] that plays no role in the present discus- 
sion. For the marginal solution q\ is fixed by equating the lowest fluctuation 
eigenvalue to zero, which gives 

l -p{p-l)p 2 J 2 qr 2 (l- qi f = l. (3) 

The condition dF/dqx = then yields x = x(T) = (p-2)(l-q 1 )/q 1 . This dy- 
namical transition sets in at temperature T c = J{p(p—2) p ~ 2 /2(p — l)^ -1 } 1 / 2 
where x comes below unity. The same transition temperature follows from 
dynamics. 



2.2 Components 

A state, called a component by Palmer, [2] is labeled by a = 1,2, •• ■ ,Af, 
and has a local magnetization profile mf = (Si) a . Its free energy F a is a 
thermodynamic potential that determines the internal energy and the entropy 
by its derivatives. In the present model F a — FTAp(mf) is know explicitly. 
It is a minimum of the "TAP" free energy functional [14] [10] [11] 

F TAp{ m i) = ~ J ii - i P m ii--- m if,- H ^2 m i 

il<---<i p i 

NT NB J 2 

— log(l - q) - -t^- (1 + (p - 1)^ - pgr- 1 ) (4) 

where q = (l/A)J]-mf is the self-overlap. The state a occurs with weight 
p a that is set by the type of experiment one describes. (In practice these 
weights are usually unknown.) Given the p a 's one can define the "component 
averages" such as F = ^ a p a F ai C = ^2 a p a C a and even the complexity 
[1] [2] X = — ^2 a p a ]nPa- For any observable, the component overage is the 
object one obtains when measuring over repeated runs and averaging over 
the outcomes. According to the Gibbs weight the probability of occurrence 
is Pa = cxp(-(3F a (T))/Z, with Z - £ a exp(-/3F a ). 

The nice thing of the present model is that many questions can be an- 
swered directly. After setting dFTAp/dmi — 0, we can use this equation to 
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express F a in terms of q a alone. This gives the simple relation F a = Nf(q a ) 
where 

Tq T 

The resulting saddle point equation for q a coincides with the marginality for 
qi given below eq. (2). Since F a only depends on the selfoverlap q a , it is self- 
averaging. In the paramagnet one has m, = q = 0, so both eqs. (4) and (5) 
reproduce the replica free energy F — —Nf}J 2 /A. From the replica analysis 
we know that at T~ the value of q a « qi is qi = q c = (p — 2)/(p — 1). The 
component free energy F a = Nf(q c ) exceeds the free energy of the param- 
agnet by an extensive amount. As expected from experimental knowledge on 
glasses, the internal energy is found to be continuous. At T~ the free energy 
difference is solely due to the lower entropy, S a = — N/3 2 J 2 /A — I c , where 

J c = 7vQlog(p-l) + ^-l) (6) 

is the value of complexity at the transition point. 

This discussion supports the picture of the glass as a disordered solid, 
where the entropy of the component the system condenses in, and thus the 
component average S, is much smaller that the entropy of the paramagnet. 
In real glasses this loss of entropy is due to the reduced phase space that 
arises by trapping of the atoms in a glassy configuration. In the quantized 
system S will vanish at T = 0. [13] 



2.3 Value of the complexity 

Kirkpatrick, Thirumalai and Wolynes [15] were the first to study the role 
of the complexity for Potts glasses in static situations in the temperature 
range T g < T < T c , see also [4] for Ising spin glasses. Statically ( that is to 
say, on timescales ~ exp(iV)) the system condenses into a state with higher 
free energy but with complexity such that the total free energy is exactly 
equal to the paramagnetic free energy. Here we will investigate the role of 
the complexity on timescales ~ iV 7 , relevant for the dynamical transition at 
T c . 

The free energies discussed for this problem are plotted in Figure 1. 

A simple calculation shows that the 'experimental' free energy of CHS and 
CS, and the marginal replica free energy obtained from eq. (2) [12] have the 
following connection with the component average free energy F = Nf(q\): 

F e xp F TX C (7) 
_ _ TT, 
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Fig. 1. Free energies of a spherical spin glass with random quartet couplings, after 
subtraction of the paramagnetic value, a) Marginal replica free energy b) Static 
replica free energy c) "Experimental" free energy, obtained by integrating C/T and 
by analysis of the degeneracy of the TAP states. 

Since x(T c ) = 1 both expressions are at T c equal to the paramagnetic free 
energy. 

In order to trace back the difference between (7) and (8) we have decided 
to redo the analysis of the TAP equations. Hereto we consider the generalized 
partition sum 

Z u = J2 e ~ Ul3Fa{T) =e~ PFu (9) 

a 

For u — we thus calculate the total number J\f of TAP-states, while for 
u = 1 we consider their partition sum. The sum over the TAP states can be 
calculated using standard approaches. [16] [14] A 1RSB pattern is assumed for 
the 6 order parameters. For instance, q a p = (l/N) mfmf takes the values 
qd for a = j3 and qi for a ^ (3 both inside a x * x diagonal block of the 1RSB 
Parisi matrix, while vanishing outside these blocks. At fixed breakpoint x the 
12-dimensional saddle point can be found explicitly. For the long time limit of 
the dynamical approach the marginality condition should be taken, [12] in the 
form given in eq. (3). As expected, the above replica expression for qi is found 
back as solution of df(qd)/dqd = at qd = qi- The result q\ = qd asserts that 
the mutual overlap between different states in the same cluster is equal to 
the selfoverlap. Like in the replica calculation of the ordinary partition sum, 
x can still take any value. In analogy with the marginal replica calculation of 
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eq (2), we expect x to be fixed by the vanishing of a fluctuation eigenvalue. 
We have therefore analyzed that 1 2 x 12 longitudinal fluctuation matrix at 
marginality. For any value of x it automatically has 3 zero eigenvalues, proving 
the marginality. Another eigenvalue vanishes for x = 0, x, = 1 (twice) and for 
x = x(T)/u. 

From this we infer that x, = x(T)/u, and thus 1 = ul c /x(T). In our case 
u = 1 it just implies that the calculated complexity is the replica value I, and 
not I c , the one of CS. This leads to the conclusion that nothing went wrong 
in the replica calculation of the dynamical phase transition: the replica free 
energy is a generating function for the mean field equations, and its saddle 
point value is the logarithm of the partition sum. 

This conclusion has been supported by a calculation for the spherical p- 
spin glass in a transverse field r. [17] In that extension of the model there 
again occurs a dynamical transition from the paramagnet to a 1RSB spin 
glass state, at transition temperature T c (r). The paramagnet of this model 
is non-trivial. There is a first order transition line (that we call pre-freezing 
line ) separating regions with large and small ordering in the z-direction. 
[18]. This line intersects the dynamical PM-1RSB transition line at a point 
(T*, r*). Beyond this point there occurs a first order PM-SG transition with 
a finite latent heat. We expect the location of the transition line to follow from 
matching of free energies. The replica free energy is indeed suited for that, 
while the 'experimental' free energy does not lead to a meaningful match. 

The free energy is F rep i is the physical one, in the sense that it takes 
into account the correct value of the complexity. Nevertheless the increase of 
complexity, I <~ 1/T for low T, remains to be explained. 



2.4 Specific heat versus energy fluctuations 

It would be nice to have a measurable quantity that probes the multitude 
of states. One object that should be accessible, at least numerically, is the 
specific heat. The standard expression C = dU/dT = ^2 a d(p a U a )/dT is likely 
to differ from the component average fluctuations of the internal energy: C = 
E a Ca = EaPadUa/dT = /3 2 E a P*( AU Z) ■ I 2 ] Thc interesting question is 
whether their difference is extensive. Based on experience in a toy model, [19] 
we think it generally is in systems with 1RSB. Since in the present model the 
energy fluctuations are too small at H = 0, [20] it can only occur in a field. 
From the internal energy in a small field we obtain 

1C(T, H) = I/? 2 J 2 (l + (p - l^-pq?- 1 ) 

_^(p-l) a <p-2)(l-«)' (10) 



On the other hand, a short calculation shows that C a = —Td 2 F a /dT 2 re- 
mains only a function of q a at the marginal point, which in the present 
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model takes the field- independent value q\ = qa- This implies that C is field- 
independent as well, thus satisfying the Parisi- Toulouse hypothesis [21]. In- 
terestingly enough, we find C < C, whereas Palmer derives the opposite at 
equilibrium. Our reversed "dynamical" inequality is a new result that is due 
to the marginality. 

The reversed dynamical inequality occurs due to non-equilibrium effects. 
We conjectured that it generally takes place outside equilibrium, for instance 
in cooling experiments above T c in the three dimensional Edwards- Anderson 
spin glass. Some numerical support for this behavior was found, see [22] 

3 Directed polymer in a correlated random potential 

We introduce a new, simple model with a static glassy transition. Consider 
a directed polymer (or an interface without overhangs) z(x) in the section 
1 < x < L and 1 < z < W of the square lattice with unit lattice constant. In 
the Restricted Solid-on-Solid approximation the interface can locally be flat 
(z(x+ 1) = z(x); no energy cost) or make a single step (z(x + 1) — z(x) = ±1; 
energy cost J); larger steps are not allowed. The polymer is subject to periodic 
boundary conditions (z(0) = z(L)) and we allow all values of z(0). 

Further there is a random energy cost V (z) per element of the polymer at 
height z. Note that this is a correlated random potential, with energy barriers 
parallel to the x-axis. 

3.1 The partiton sum 

The partition sum of this system can be expressed in the eigenvalues of the 
tridiagonal transfer matrix T that has diagonal elements exp(— (3V(z)) and 
off-diagonal elements exp(— [3J) 

w 

Z = tre-W = trT L = ^ {A W ) L (11) 

w — 1 

For a pure system (V(z) = for all z) at temperature T = 1/(3 Fourier 
analysis tells that for small momentum A(k) = 5?(1 + e -P J + ik y (i - e -pj+ik^ 
« exp[— [5Jb — Dk 2 / (2tt 2 )] with bulk free energy density Jb and diffusion 
coefficient D that can be simply read off and are temperature dependent. 

We shall consider the situation of randomly located potential barriers 
parallel to the x-axis. Hereto we assume binary disorder, so V(z) = with 
probability p = cxp(— /i) or V(z) = V\ > with probability 1 — p. Eq. (11) 
is dominated by the largest eigenvalues. It is well known that they occur due 
to Lifshitz-Griffiths singlaritics. These are due to lanes of width t 3> 1 in 
which all V(z) = 0, bordered by regions with V(z) ^= 0. These dominant 
configurations are the "components", "TAP states" or "metastable states" 
of our previous discussion. The eigenfunction centered around z has inside 



Complexity as the driving force for glassy transitions 



9 



the lane the approximate form cos [71- (z — z )/£] while it decays essentially 
exponentially outside due to the disorder. These states can thus be labeled 
by a = (z a ,£ a ). Since k — > ir/£ the free energy of this state follows as 

PF e = -LlnA e ^Pf B L+^ (12) 

The typical number of regions with £ successive sites with V = is Me = 
W(l —p) 2 p e . We now choose W — exp(AL 1//3 ) so the states with width £ have 
a configurational entropy or complexity Te = \nMe « AL 1 / 3 — fj,£. 



3.2 The TAP-partition sum 

For large L we may restrict the partition sum to these dominant states. We 
thus evaluate, instead of eq. (11), the 'TAP' partition sum 

Z = Y t *te- f,Fl (13) 

1 

Note that it is obtained by simply omitting the contributions of states with 
low eigenvalue (high free energy). The total free energy 

0F = - In Z = pf B L - XL 1 / 3 + fi£ + ^ (14) 

has to be optimized in £. The largest £ which occurs in the system can be 
estimated by setting Me ~ 1, yielding 

XL 1 / 3 

^max (15) 

It is a geometrical length, independent of T. Let us introduce D = Dp? /A 3 . 
The free energy of this state reads 

(3F = Pf B L+^XL l l 3 D (16) 

At low enough T the optimal length is smaller than £ ma x, 

1/3 



-if) 



(17) 



The free energy of this phase is 



(3F = (3f B L + iAL 1 / 3 (3^ 1 / 3 - 2) (18) 

For D > 1 (T > T g ) the interface is in an essentially non-degenerate state. 
For D < 1 it lies in one of the Me 3> 1 relevant states, which is reminiscent 
to a glass. So the model has a glassy transition at D = 1. 
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The internal energy of a state of width I is 




AL\ 1/3 18D 



(19) 



At D = 1 this coincides with the paramagnetic value, simply because I — > 
(■max- It is easily checked that the free energy (18) is a thermodynamic po- 
tential, and yields the same value for U . At the transition it branches off 
quadratically from (16). In the glassy phase the specific heat 



exceeds the component averaged specific heat C = Lcb + (£/2^ 2 )<9t<9/3-D. In 
contrast to previous model, the specific heat it is larger in the glassy phase 
than in the paramagnct. This is because the free energy is lower. 

3.3 On overlaps and hierarchy of phase space 

In a given realization of disorder we define the 'overlap' of two states a and 
b, centered around z a and Zb, respectively, as 



In the high temperature phase there is one non-degenerate state, so P(q) = 
S(q — qi). In the glassy phase we expect that q a b = Qi for all optimal states 
(a, b) at temperature T . The reason is that at thermodynamic equilibrium 
the whole phase space can be traversed, and negligible time is spent in non- 
optimal states. If so, then though there are many states, one still has P{q) = 
S(q — qi) and there is no replica symmetry breaking. This is standard for 
equilibrium situations without frustration. 

This puts forward the picture of replica symmetry breaking and hierarchy 
of phase space being a dynamical effect. At given timescale only some nearby 
states can be reached, "states within the same cluster" . At larger times other 
clusters can be reached, and for times larger than the ergodic time of a large 
but finite system, all states are within reach. Only in the thermodynamic 
limit phase space splits up in truely disjoint sets. To investigate the validity 
of this picture in detail, one should solve the dynamics of the polymer. 

3.4 The polymer model at T = 1/T 

The comparison to the p-spin model is most direct when we compare the 
p-spin model at temperature T with the polymer model at temperature 
T = 1/T. In this interpretation, coming from high T, the polymer under- 
goes a gradual freezing into TAP states. This truely becomes relevant when 
the domain size is of order £ ~ L 1 / 3 , where the complexity starts to be smaller 



° = % = CBL + W° TdpD + r 3 ( XL D- 5 ) 1/3 (Td T D) 2 (20) 



Qab = Hm {8 z (p), z J z (t),z b ) 



(21) 
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than log IT = XL 1 / 3 . This gradual freezing shows explicitly that the dynam- 
ical transition, as found in the mean field p-spin glass, is smeared in finite 
dimensions. 

For T going down to T c , the polymer gets captured in states with free 
energy closer and closer to the lowest free energy state available at that 
temperature. 

As in 1RSB spin glasses, the complexity also vanishes to leading order for 
T I T c . In the low T phase the complexity is no longer of order L 1 / 3 . This is 
similar to the low T phase of the static p-spin model, where the complexity 
is non-extensive in the glassy phase. 

When considered as function of T, the specific heat makes a downward 
jump when cooling the system from large T below T c . The absence of a sharp 
dynamical transition and the vanishing of the complexity that occurs in this 
polymer model as function of T are very analoguous to the expected behavior 
of realistic glasses. 
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